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Abstract 

We discuss experimental aspects related to the e + e~ — > ZZ process and to 
the search for anomalous ZZV couplings (V= Z, 7) at LEP2 and future e + e~ 
colliders. We present two possible approaches for a realistic study of the 
reaction and discuss the differences between them. We find that the optimal 
method to study double Z resonant production and to quantify the presence 
of anomalous couplings requires the use of a complete four-fermion final-state 
calculation. 
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INTRODUCTION 



Pair production of Z bosons is one of the new physics processes to be studied at LEP2 
and at future high energy e + e~ colliders. Although it is a process with a rather low cross 
section (below 1 pb) and experimentally difficult to observe (large and almost irreducible 
backgrounds), LEP2 gives the first opportunity to perform a measurement and to look for 
deviations from the Standard Model (SM). In addition, a good understanding of the process 
is necessary, since it is one of the relevant backgrounds in the search for the Higgs particle. 
At future e + e~ colliders, with luminosities of the order of 100 fb _1 , several thousands of 
events will provide stringent tests of the SM. 

The study of triple-gauge boson couplings is one of the key issues at present and future 
colliders. Anomalous Z7V couplings have been searched for at the Tevatron and at LEP 
0]. The first experimental limits on anomalous ZZZ couplings have been provided by the 
L3 Collaboration @. 

The report is organized as follows. First, the SM amplitude is presented and the effects 
of possible anomalous ZZV couplings at LEP2 and at future e + e~ colliders are discussed. 
Second, we describe two reweighting approaches developed for the search for anomalous 
couplings at LEP2. These approaches are compared and their differences are pointed out. 
Finally, the fitting techniques used to quantify the possible presence of anomalous couplings 
are briefly presented. 

STANDARD MODEL AMPLITUDE FOR THE e+e~ ZZ PROCESS 

The diagrams contributing at first order to the e + e _ — > ZZ process in the Standard 
Model are shown in Figure |]. We will assume a collision in the center-of-mass system with 
total energy yfs and neglect the effect of the electron mass. The following notation is used: 

• Electron four-momentum k and helicity a: 



a 6 {-1,1} 



Positron four-momentum k and helicity a: 



t. — (VI _Vl £ 

a I 2 ' 2 Cz 



<f €{-1,1} 



Z four-momentum g Zl and polarization e Zl : 



g Zl = (E, y/E» - Ml g); 
e Zi = ezi(AzJ; 



A Zl e {-1,0,1} 



I z l z 2 



Z four-momentum gz 2 and polarization ez 2 : 



Tpl v I z 2 L l 



9z 2 = (•£"> — y 2 — M| 2 g); 
ez 2 = ez 2 (Az 2 ); Az 2 G {-1,0,1} 

where the electron is assumed to collide along the +z axis (e 2 ), and the Z -with mass M Zl - 
goes along the direction given by g = (sin#z cos^z, sin^z sin^z, cos#z)- The masses M Zl 
and Mz 2 are not assumed to be equal to the on-shell mass m z because in the following they 
will be consider as virtual particles decaying into fermions. 

The matrix element for the e + e~ — > ZZ reaction is determined by the same method 
followed in III. It reads: 



M zz (a,cr, A Zl , A Z2 ) = -{g„ 



Zc+c~ \2 



S(e* Z2 ,q Zl ,e* Zl ,a) S{e Zl ,qz 2 ,e* Z2 ,a) 



-2(kq Zl ) + Ml -2(kq Z2 ) + Ml 



The functions S(e a , q b , e b , a) are given by: 

x ( 

S{e a , q b , e b , +) = ( £ + ^ - e ° - e 3 a ' 



" -9? -9b> -Qb+nl ^ 

\ -q\ - ~ib + 1! J 



F° - f 3 



-ei - let 



S(e a ,q b ,e b ,-) = e ° a + e 3 a , e\-%e\ 



-96 + Qb, Qb- Hb 

A _l „-„2 rz „o „3 



\ Qb+Wb, V« - 96 - Qb J 



e b e b 



(2) 



(3) 
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where the components of the four- vectors are denoted by superscripts. The left/right effec- 
tive couplings of fermions to neutral gauge bosons are given by: 



gf = -2Q f sin^w (^2G>|) (4) 

9 ^ = gf + 2h(V2G,ml) 1/2 (5) 

gf = Q f (4vra(M 7 2 ,)) 1/2 (6) 

gf = gf (7) 



where Qf is the charge of the fermion f in units of the charge of the positron, and the 
electromagnetic coupling constant a(M 2 *) is evaluated at the scale of the virtual photon 
mass M^». I 3 is the third component of the weak isospin (±1/2), sin 2 #w is the effective 
value of the square of the sine of the Weinberg angle and is the value of the Fermi 
coupling constant. The effective couplings to the Z absorb the relevant electroweak radiative 
corrections at the scale of the Z [|]. They are obtained by the substitutions: 

sin 2 w — > sin 2 #w (8) 

a ■ 2/ 2 n - V2G,ml (9) 
4 sin &w cos^w 

The experimental signature of a e + e~ — > ZZ process is a final state with four fermions, 
due to the unstability of the Z particle. A distinctive feature is that the invariant masses 
of the two pairs, ff and f'f , peak at the Z mass mz- The angular distribution of the decay 
products keeps information on the average polarization of the Z boson. In addition, the Z 
decay amplitude has to be included for a correct treatment of spin correlations. Assuming 
that fermion masses are negligible compared with mz, this amplitude is given in the rest 
frame of the Z by: 

M Zif{ (\ Zi ,X,\)=g z x ff M Zi 5 x> _j [e Zi {v x -i\v 2 )\ (10) 



v\ = (0, cos Of cos <j)f, cos9f sin 0/, — sin Of) 
V2 = (0, — sin 0/, cos0/, 0) 



(11) 
(12) 



where p = (sin0/ cos0/, sin 9f sin 0j, cos8f) is the direction of the fermion momentum and 
A, A are the helicities of fermion and antifermion, respectively. 

ANOMALOUS COUPLINGS IN THE e+e~ -> ZZ PROCESS 

Anomalous gauge boson couplings lead to interactions of the type shown in Figure |2|. 
The coupling ZZV, with V = Z or 7, does not exist in the Standard Model at tree level. 
Only two anomalous couplings are possible if the Z bosons in the final state are on-shell, 
due to Bose-Einstein symmetry. In principle, five more couplings should be considered if at 
least one of the Z bosons is off-shell. However, as discussed in Appendix [A], the new terms 
must be of higher dimensionality and are suppressed by orders of m z r z /A 2 , where A denotes 
a scale related to new physics. We will concentrate on the most general expression of the 
anomalous vertex function at lowest order 131: 



r^ v = { ifX [fazx + qz 2 ) a 9^ + (g Zl + qz 2 r 9^] + if? ^ P (fe - qz 2 ) P } (13) 



mi 

A non-zero value of fj leads to a C- violating, CP- violating process, while terms associated 
to /g are P-violating, CP-conserving. Using again the formalism followed in || we obtain 
the explicit expressions for the anomalous contributions: 

M f / C (a,a, A Zl , A Z J = - i etfg^ ^ [e£(e£ + ioe%) + ej^eg + iaefj] (14) 
M f / C (a, a, A Zl , A Za ) = - i egg^ ^ 8^ (e 1 ^ + ^ e* Zia e* Za/3 (q Zlp - q z . 2P ) (15) 

Note that no (s — my) factors are present in the final expressions. Compared to the SM 
amplitude all anomalous contributions increase with the center-of-mass energy of the colli- 
sion. We want to bring the attention to the fact that these anomalous ZZ7 couplings are 
different from those considered in the e + e~ — > Z 7 anomalous process ||. Therefore, not 
only the two anomalous ZZZ couplings, but all four anomalous parameters remain essentially 
unconstrained at present. 

Anomalous ZZV couplings manifest in three ways: 
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• A change in the observed total cross section e + e~ — > ZZ. 

• A modification of the angular distribution of the Z. 

• A change in the average polarization of the Z bosons. 

The effect of anomalous couplings in the e + e~ — > ZZ process at Born level is illustrated in 
Figures ^|and [|for the center-of-mass energies of \/s = 190 GeV and y/s = 500 GeV, respec- 
tively. The anomalous distributions are determined for the values — 3; % — 4, 5; V — Z, 7. 
Both CP-violating and P-violating couplings are found to produce a global enhancement in 
the number of events. This increase is very clear at \/s = 500 GeV. There are moderate 
changes in the angular shape at y/s = 190 GeV. At \/s = 500 GeV the situation is different. 
The copious anomalous production at large polar angles starts to compensate the huge peaks 
of the SM differential cross section at low polar angles. The SM divergent behaviour happens 
in the limit mz <C \/s, where the process tends to a t-channel process with production of 
massless bosons. 

Anomalous couplings also modify the average polarizations of the Z bosons, as shown in 
Figures [?] and [7[ At y/s = 190 GeV the observed change depends on the particular size 
and type of the anomalous coupling under consideration. For CP- violating couplings there is 
always an increase of the production of bosons with different polarizations (longitudinal ver- 
sus transverse). At \fs = 500 GeV all couplings show a similar behaviour: an enhancement 
of longitudinal-transverse production and a suppression of transverse-transverse production. 
This is an interesting feature, since the SM process has the opposite behaviour. The fraction 
of states in which the two bosons are longitudinally polarized is below 0.5% at these energies, 
both for SM and for anomalous production. What is physically observable is a modification 
of the angular distributions in the center-of-mass frame of the Z — ► ff decays: in absence of 
anomalous couplings, both Z decays will proceed preferentially along the direction of the Z 
momenta, whereas one of the decays will preferentially occur at 90° if the process is highly 
anomalous. 
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Summarizing, at energies close to the threshold of ZZ production the sensitivity to 
anomalous couplings is weak. For an integrated luminosity of 200 pb -1 , tens of events 
are expected to be selected. The main anomalous effect is an increase in the cross section, 
and one expects small improvements from the variations in the angular distributions. At 
higher energies, with luminosities of the order of 100 fb _1 , all anomalous effects contribute 
coherently to enhance the sensitivity: a huge increase in the cross section, especially at 
large polar angles, and a clear correlation between the angular distributions of the Z decay 
products. 

INCLUSION OF ANOMALOUS COUPLINGS. REWEIGHTING PROCEDURE 

In order to take into account anomalous effects with sufficient accuracy, the correct 
matrix element structure has to be implemented. In many cases, and from the practical 
point of view, the generation of events for different values of anomalous couplings is not 
convenient. A more attractive method is to set up a procedure to obtain the Monte Carlo 
anomalous distributions as a function of a single set of generated events. This is the role of 
reweighting methods. 

Let us consider a set of events generated according to the Standard Model differential 
cross section. New distributions taking into account the anomalous couplings are obtained 
when every event is reweighted by the factor: 



2 

J2 (M zz + M A c) M Zlff M Zaf , f 

Az x ,Az 2 



2 

]T M zz M Zlff M Z2f ,p 

Azj^ >Az 2 



(16) 



The weight W zz (a, A, A'; 0) depends on the helicities of the initial electron (er) and of 
the final fermions (A, A'). It also depends on the kinematic variables defining the phase space 
(Q). For convenience we choose the following set: 

• The invariant masses of the f f and f f ' systems: M Zl , M Za . 



The polar and azimuthal angles of the ff system : 8 Z , <pz- 



• The polar and azimuthal angles of the fermion f after a Lorentz boost to the rest frame 
of the ff system: 9f, 0/. 

• The polar and azimuthal angles of the fermion f after a Lorentz boost to the rest 
frame of the ff system: 9f, <f)f. 

The previous result can be extended to take into account other non-resonant diagrams 
like e + e~ — > Z7* — > fff f and e + e~ — > 7*7* — > ffff . These diagrams can not be neglected 
for a correct analysis of double Z resonant production ||. The final weight is: 

V Zl (Ml)V Z2 (Mi 2 ) M ac M Zlf{ M Z2f , f , 



£ V Yo (Ml)V v ,(Ml) M VoV a M Vofi M v , f , f , 



V Q ,V' 



(17) 



where a sum on all intermediate V G {Z,7*} is assumed. The propagator factors T>v are 
defined as follows: 

1 



q 2 — m| +« T z q 2 /m z 



V,(q 2 ) = 3 



(18) 
(19) 



where the imaginary component takes into account the energy dependence of the Z width 
around the resonance. The expressions for M Zl *, M 7 * 7 * , M » f j are obtained by the same 
method used for M zz and M Zf j. Explicitly, they can be obtained by substituting Z by 7* 
where necessary: 



ez(A z ) -> e r (A 7 ) 



M Zl 



,Zff 



5+ 



Tff 

9- 



(20) 
(21) 
(22) 
(23) 



Weights according to W YV (a, X,X';Q) have been implemented in a FORTRAN program. 
The approach is well suited for events generated with the PYTHIA e + e~ — > Z/7* Z/7* — > 



fff'f' generator |J. This implementation will be identified as the "NC08 approach", because 
all eight neutral conversion diagrams are considered. Several checks have been done in order 
to make sure that the calculations are correct. There is also agreement with the results 
obtained in and in ||. 

After reweighting, distributions according to given values of the anomalous couplings 
fii /jT are obtained. Detector effects are correctly taken into account if events are reweighted 
at generator level. 



INITIAL STATE RADIATION EFFECTS 

There are several references providing valuable information on the e + e~ — ► ZZ process. 
A specific SM generator for e + e~ — > Z/7* Z/7* (7) — > fff'f '(7) without anomalous couplings 
exists in PYTHIA 0. The calculation reported in the previous section is well suited for this 
MC generator, but initial state radiation effects (ISR) need to be taken into account. We 
assume that the differential cross section can be expressed as follows: 

da(s) r , da(s') 



I ( ' ) rffPhase Space') ^ ^ 



<i(Phase Space) J <i(Phase Space') 

where cr(s') is the (undressed) cross section evaluated at a scale s', and a(s) is the cross 
section after inclusion of ISR effects. The radiator factor R(s, s') is a "universal" radiator, 
that is, independent of specific details of the matrix element. With this assumption ISR 
effects are accounted for by evaluating the matrix element in the center-of-mass system of 
the four fermions, at the corresponding scale s'. 

In Reference a specific e + e~ — > ZZ(7) — > fff'f'(7) generator for anomalous couplings 
studies is presented. It takes into account ISR effects with the YFS approach up to 0(a 2 ) 
leading-log. It has some limitations, like the absence of conversion diagrams mediated by 
virtual photons. 

The Standard Model cross section (f¥ — ftf = 0) from ]7| shows agreement at the 
percent level with the one determined in Reference |J, where it is shown that all significant 
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radiation effects come from "universal" radiator factors. This implies that an approach 
based on Equation El is justified in terms of the required precision. 



THE COMPLETE e+e~ fff'F PROCESS 

Additional non-resonant diagrams are taken into account in SM programs for general 



four-fermion production, like EXCALIBUR [10|. Under a reasonable set of kinematic cuts, 



the relative influence of those diagrams can be reduced, but not totally suppressed. This is 
due to the low cross section for resonant ZZ production. Typical examples are those involving 
charged currents (relevant in e + e~ — > UqUJI, e + e~ — > uddu, . . . ) or multiperipheral effects in 
e + e~ — > e + e~ff . In addition, the influence of Fermi correlations in final states with identical 
fermions (e + e~ — > ffff) is unclear. 

In order to include these effects, the EXCALIBUR program has been extended. All 
matrix elements from conversion diagrams with two virtual Z particles M zz (a, A, A'; Q) are 
modified in the following way: 

]T M AC M Zlfi M Z2t ,p 
A zz (er, A, A'; fl) = ^ (25) 

M| z xc (a, A, A'; Cl) -> M| z xc (a, A, A'; 0) (l + A zz (a, A, A'; 0)) (26) 



where Mzz, Mac, M Zli j and M Z2i ,p are the same terms defined in the NC08 approach. Based 
on this modification it is straightforward to define an alternative reweighting. It will be 
identified as the "FULL approach" in the following. More detailed studies are reported in 
the next section. 



THE NC08 APPROACH VERSUS A FULL TREATMENT 

All checks presented in this section require a precise definition of a "ZZ signal". Channels 
involving electron or electronic neutrino pairs in the final state have a non-negligible contri- 
bution from non-conversion diagrams. Also final states with fermions from the same isospin 
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doublet ((£, i^),(u,d),(c,s)) show a non-negligible charged current contribution. Therefore 
stringent cuts must be applied in order to select a sensible e + e~ — > ZZ experimental signal. 
Our signal definition implies the following cuts: 

• The invariant masses of the two fermion-antifermion candidate pairs must be in the 
range 70 GeV-105 GeV. 

• In the final states with electrons, these electrons must verify | cos^l < 0.95. 

• In the final states with WW contributions, the invariant masses of the fermion pairs 
susceptible to come from W decay must be outside the range 75 GeV-85 GeV. 

The SM cross sections within signal definition cuts at \/s = 190 GeV are shown in Figure |8| 
for the different four-fermion channels. The EXCALIBUR generator is used. Two determi- 
nations are shown: one taking into account all the Standard Model diagrams and the other 
considering the neutral conversion diagrams only. Note that, in some cases, there are large 
differences between the two calculations. This already points to the convenience of using a 
full four-fermion approach, even in the presence of strong cuts. 

The next study is devoted to the SM matrix elements. For the same set of neutral 
conversion diagrams, the results obtained by EXCALIBUR and by the NC08 approach are 
compared. The relative differences are shown in Figure |9|. Two groups are considered. The 
first group corresponds to all processes without Fermi correlations, that is, those in which 
there are no identical fermions in the final state. A perfect agreement is observed in this 
case. The second group contains those processes in which there are identical particles in the 
final state. Let us note that no effort has been done in the NC08 approach to antisymmetrize 
the matrix element in the presence of identical fermions. Although the effect is understood 
and it can be trivially included, we want to show that it is not totally negligible. It may 
exceed 10% for some phase space configurations. 

Figure [II] shows the distribution of weights at a/s = 190 GeV for non-zero values of the 
anomalous couplings. The fact that the distributions are not extremely narrow indicates 
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the presence of effects other than just an excess of events. To compare the implementations 
in the presence of anomalous couplings, the Standard Model distributions are reweighted 
according to the NC08 and FULL approaches. Again, only neutral conversion diagrams are 
considered. The relative differences between the weights assigned by the two approaches are 
shown in Figure |TT|. There is agreement at the percent level. Fermi correlations in the final 
state have a small effect on the weights. The reason could be related to the fact that the 
discrepancy factorizes in a similar way SM and anomalous terms. 

The last study evaluates the influence of non-conversion diagrams in the presence of 
anomalous couplings. The averages of the reweighting factors for the FULL and the NC08 
approaches are compared in tables |- |I[ The set of cuts defining the ZZ resonant region 
is applied in all cases. The differences are typically below 10%, but not negligible. This 
points again to the convenience of considering all possible diagrams contributing to the 
e + e~ — > fff'P process. 



MEASUREMENT OF ZZZ AND ZZ7 ANOMALOUS COUPLINGS 

To determine the values of the anomalous couplings from the data, the histogram of the 
most relevant variable for each four-fermion channel may be used. The following binned 
likelihood function is then maximized: 

log(£) = £ [N data {j) log N expected (j; fl) ~ N expected (j; fl)} (27) 
i=i 

The expected number of events is computed as: N expected (j) fj) = N signa i(j; fl) + 
Nbackground(j) ■ The background contribution does not depend on the anomalous couplings. 
The signal contribution includes all four-fermion final states compatible with the exchange 
of two Z bosons. It is computed by reweighting the Standard Model distributions with the 
FULL approach and taking into account the fY values of the anomalous couplings. 

This method has been applied in the analysis of e + e~ — > fff'P final states by the L3 
Collaboration H. The discriminating variables are the invariant masses of the lepton pairs 



12 



in leptonic decays and neural net outputs for hadronic decays. They obtain the following 
95% confidence level limits on the existence of anomalous ZZV couplings: 

-1-9 < ff < 1.9; -5.0 < fi < 4.5; -1.1 < ft < 1.2; -3.0 < f 5 < < 2.9. 

As an example of the use of more sensitive variables at higher energies we will consider 
the semileptonic process e + e~ — > £ + £~qq at y/s = 500 GeV. This channel is expected to 
provide a clean signature for ZZ production if the Higgs mass is away from the Z mass 
region. In order to enhance the ZZ signal component, the invariant mass of the leptons is 
required to be in the 70 GeV to 150 GeV range, the recoiling hadronic mass has to be larger 
than 50 GeV and the polar angle of electrons and positrons, 9 C , must satisfy | cos# e |< 0.8. 
The differential cross section of the process in the presence of an anomalous coupling / can 
be expressed as follows: 



d 2 a 



d(Oi) d(0 2 ) 



d 2 a 



d(Oi) d(0 2 ) 



(l + fOr + f 2 ) (28) 

/=o 



The variables 0\ and 2 are functions of the phase space variables of an event. They 
are independent of /. The previous equation guarantees that the maximal information on 
/ is obtained by a study of the event density as a function of the variables 0\ and 2 . 
These variables are usually called "optimal observables" . Given an event characterized by 
the phase space point 0, a simple expression for the optimal observables is: 

om= gmif ± mi = -!)_! (30, 

2 

where W (0; /) is the reweighting factor used to transform SM distributions into anomalous 
distributions at Cl for an anomalous coupling /. At ^J~s = 500 GeV and for small values 
of / the dominant term in the ZZ differential cross section is the one associated to 0\. In 
real life, the exact values of the phase space variables are unknown, but we may use as an 
approximation the phase space variables reconstructed by the detector. For this exercise 
we assume energy resolutions of 5% for quarks and leptons and a jet angular resolution 
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of 40 mrad. The weight is symmetrized under the interchange of quark types, assuming 
that the quark flavour can not be identified. The SM distribution of the variable 0\ for 
an integrated luminosity of 100 fb^ 1 is shown in Figure O, together with the effect of an 



anomalous coupling /f = 0.01. Note that the ratio between the anomalous and the SM 
distributions is not constant. This is a demonstration of the sensitivity of the variable and 
of the existence of anomalous effects different from a simple change in the total cross section. 
From a likelihood fit to the anomalous distribution we obtain a value of f§ = 0.010 ± 0.002. 
Note the increase in sensitivity as compared to LEP2 present limits. 
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APPENDIX A: MOST GENERAL ANOMALOUS ZZV COUPLINGS 

The most general Lorentz invariant anomalous coupling structure for the ZZV vertex 
function is, similarly to the WWV case 0: 

r z?zW = { ifl ((fe + lz 2 ) a g* + (g Zl + gz 2 y <T) + ifl - qz 2 ) P } 

+ 



•z - - " - - - - - - ' ■■■ 



+ fl ((? Zl + qz 2 ) a - (q Zl + qz 2 f gn - fl e al3w (? Zl + qz 2 ) P 

- $r (gz 1 ~ qz 2 T ^ pu (qz 1 + qz 2 ) P (g Zl - qzX } (Al) 
The global factor (s — my)/m| is introduced by convention [] in order to preserve gauge 



1 The presence of mz in the denominator is arbitrary. It allows the introduction of a dimensionless 
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invar iance when V = 7 and Bose-Einstein symmetry for V = Z in the on-shell limit. Note 
that only the terms associated to f[ and fY survive in the limit in which both Z are on-shell. 

The ZZV vertex function must be symmetric under the interchange (Z\,a) <-> (Z2,/3). 
For the terms associated to fY, f%, f%, f% and fj this requirement forces the presence of an 
additional Lorentz invariant factor: any antisymmetric function of (q z ~Qz 2 )- Our ininiixial 
choice is: (q z — g| 2 )/A 2 , where A is a scale related to new physics. Two comments are in 
order here: 



The anomalous couplings f{ , / 2 , / 3 , / 6 , fj are necessarily associated to Lagrangians 



V fV 



of higher dimension than those associated to / 4 , f 5 



• The sensitivity to the anomalous /i 7 , /J 7 , Z^, /(T? fY couplings is further reduced due 
to the relatively small size of the Z width: (q Zi — g| 2 )/A 2 ps C(rzmz/A 2 ). 

The existence of additional (qf — m|) factors when at least one of the final Z is off-shell 
has also been noticed in Reference || . For completeness, we provide the full expressions for 

fV fV 

all possible anomalous couplings matrix elements, M^ c = MJ c (a,a, X Zl , Az 2 ): 



f v 

fV 



M J3 
m AC 



M J A 6 C 



1 efr 9 Y ee 8*,-* 



,2* - 



2* * 



i efY 9 Y ee £ V* (e 1 ^ + ioeW) e* Zia (q Zlf) - q z 



'z 



2pJ 



V„Vee s 3 / 2 (gj-gSJ 



A 2 



0* c 0* 

-z 2 



fV Vee s ( q Z\ ^Zj) r 

--V n Vee s (gZj/^laj 



2 



e Zl (4 2 +^ e z 2 



4(4 + ^4 



/■V Vee s ^Zi ^z 2 / r _ / la/30 , 7 w2a/30N * * 



M^c = e 



fV n Vee s (ggi g Z 2 ) c 
J7^ ^ A 2 ' ^ 



7z a -?z a ) + «K?Z a -?z 2 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(A7) 



e ^o P (q Zl -q Z2 ) p e Zia e* Z2 p (A8) 



coupling constant without adding new unknown scale parameters. A physically more motivated 
choice is to substitute m z by a scale of new physics A. In this way, the unknown coupling is of the 
type //A 2 , with / of order unity, but the higher dimensionality of the term is exhibited. 
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TABLES 







/f = l 


Final state 


NC08 




FULL 


NC08 




FULL 


qqqq 


1.128 




1.123 


1.023 




1.022 


qqyv 


1.131 




1.150 


1.019 




1.024 


qqlT 


1.105 




1.091 


1.013 




1.006 




1.101 




1.082 


1.025 




1.013 


ml 


1.128 




1.054 


1.024 




1.007 


VVVV 


1.084 




1.170 


1.022 




1.036 



TABLE I. Average reweighting factors in the presence of anomalous ZZZ couplings. The case 
with only neutral conversion diagrams (NC08) is compared to the complete approach with all 
Feynman diagrams (FULL). 





/7 = 1 


/ 5 7 = 1 


Final state 


NC08 




FULL 


NC08 




FULL 


qqqq 


1.350 




1.336 


1.041 




1.041 


qqyv 


1.339 




1.385 


1.047 




1.050 


qqlT 


1.273 




1.235 


1.032 




1.026 


\\vv 


1.218 




1.181 


1.031 




1.004 


fflT 


1.380 




1.161 


1.033 




1.019 


vvvv 


1.254 




1.463 


1.071 




1.082 



TABLE II. Average reweighting factors in the presence of anomalous ZZ7 couplings. The 
case with only neutral conversion diagrams (NC08) is compared to the complete approach with all 
Feynman diagrams (FULL). 
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FIGURES 




FIG. 2. Diagram with anomalous ZZ7 and ZZZ couplings contributing to the e + e — ► ZZ 
process. 
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FIG. 3. Effect of non-standard couplings in the e + e~ — > ZZ process at y/s = 190 GeV. A 
collision in the e + e~ center-of-mass system is assumed. The angle Oz is the polar angle of one of 
the Z bosons and da j d cos 9 z is the differential cross section. 
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FIG. 4. Effect of non-standard couplings in the e + e~ — > ZZ process at y/s = 500 GeV. A 
collision in the e + e~ center-of-mass system is assumed. The angle Oz is the polar angle of one of 
the Z bosons and da j d cos 9 z is the differential cross section. 
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FIG. 5. Proportion of events in which both Z are longitudinally polarized (top left), both Z 
are polarized transversely (top right) or one is longitudinally and the other is transversely polarized 
(bottom) in a e + e~ — > ZZ process at = 190 GeV. Longitudinal and transverse polarizations 
are denoted by the labels "L" and "T", respectively. Predictions for SM (continuous line) and 
anomalous ZZZ couplings (dashed lines) are shown. A collision in the e + e~ center-of-mass system 
is assumed. The angle 6z is the polar angle of one of the Z bosons. 
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FIG. 6. Proportion of events in which both Z are longitudinally polarized (top left), both Z 
are polarized transversely (top right) or one is longitudinally and the other is transversely polarized 
(bottom) in a e + e~ — > ZZ process at = 190 GeV. Longitudinal and transverse polarizations 
are denoted by the labels "L" and "T", respectively. Predictions for SM (continuous line) and 
anomalous ZZ7 couplings (dashed lines) are shown. A collision in the e + e~ center-of-mass system 
is assumed. The angle 6z is the polar angle of one of the Z bosons. 
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FIG. 7. Proportion of events in which both Z are transversally polarized (left), or one is longi- 
tudinally and the other is transversely polarized (right) in a e + e~ — > ZZ process at \/i = 500 GeV. 
Longitudinal and transverse polarizations are denoted by the labels "L" and "T" , respectively. The 
two figures at the top (bottom) illustrate the effect of possible anomalous ZZZ (ZZ7) couplings. 
The Standard Model predictions are shown with a continuous line, whereas the anomalous predic- 
tions are shown with dashed lines. A collision in the e + e~ center-of-mass system is assumed. The 
angle 9z is the polar angle of one of the Z bosons. The case in which both Z are longitudinally 
polarized is not shown. It accounts for less than 0.5% of the total. 
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FIG. 8. Cross sections at yfs = 190 GeV for the different four-fermion channels taking into 
account all the diagrams (histogram) and only the neutral conversion ones (NC08). Some cuts 
(described in the text) have been applied in order to enhance the ZZ resonant contribution. Other 
diagrams are important when electrons, electronic neutrinos or fermions from the same isospin 
doublet are present in the final state. 
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FIG. 9. Distributions of the relative difference, R, between the SM squared matrix elements 
computed by EXCALIBUR and by the NC08 approach. Only neutral conversion diagrams are 
taken into account. The plot a) shows the case of final states without identical fermions (fff'f ). 
is shown in the plot. The plot b) contains those final states with identical fermions in the final 
state (ffff). The effect of Fermi correlations (no effort has been done to include them in the NC08 
approach) is visible as a tail in the distribution. 
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FIG. 10. Distributions of the reweighting factor, W, in the e + e~ — > ZZ(-y) — > fff'f'(7) process 
for different values of anomalous ZZV couplings. The reweighting factors are obtained with the 
FULL approach at y/s = 190 GeV. Only neutral conversion diagrams are taking into account. 
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FIG. 11. Distributions of the relative difference, R, between the weights computed with the 
FULL and the NC08 approaches. The process is e+e" -► ZZ( 7 ) -► fff f (7) at = 190 GeV for 
different values of anomalous ZZV couplings. Only neutral conversion diagrams are taking into 
account. 
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FIG. 12. Distributions of the optimal variable 0\ in the Standard Model (histogram) and in 
the presence of an anomalous coupling /f = 0.01 (points). The process under consideration is 
e + e~ — ► ^ + £~qq at yfs = 500 GeV and for an integrated luminosity of 100 fb _1 . The second 
plot shows the ratio between the anomalous and the SM distributions. A fit to the anomalous 
distribution leads to the value / 5 Z = 0.010 ± 0.002. 
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